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Abstract 



Initial data for solutions of Einstein's gravitational field equations cannot 
be chosen freely: the data must satisfy the four Einstein constraint equations, 
^"v We first discuss the geometric origins of the Einstein constraints and the role 

the constraint equations play in generating solutions of the full system. We 
then discuss various ways of obtaining solutions of the Einstein constraint 
equations, and the nature of the space of solutions. 
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Yvonne Choquet-Bruhat's epic work of over 50 years ago shows that if a set of 
smooth initial data which satisfies the Einstein constraint equations is given, then 
we can always find a spacetime solution of the Einstein equations which contains 
^jy an embedded hypersurface whose metric and second fundamental form agree with 

the chosen data. In the years since then, arguably the most important method 
for constructing and studying solutions of Einstein's equations has been the initial 
value (or Cauchy) formulation of the theory, which is based on Yvonne's result. 
Especially now, with intense efforts underway to model astrophysical events which 
produce detectable gravitational radiation, the Cauchy formulation and numerical 
efforts to implement it are of major interest to gravitational physicists. 

To understand the Cauchy formulation of Einstein's theory of gravitation, we 
need to understand the constraint equations. Not only do the constraints restrict 
the allowable choices of initial data for solutions; they also effectively determine the 
function space of maximally globally hyperbolic solutions of the theory, and they 
play a role in generating the evolution of the initial data via their appearance in 
the Hamiltonian for Einstein's theory. 

The goal of this review paper is to provide some measure of understanding of 
the Einstein constraints. We start in Section 2 by explaining the geometric origin 
of the constraint equations. To do this, we discuss 3 + 1 foliations of spacetimes, the 
Gauss-Codazzi decompositions of the curvature, and the consequent 3 + 1 projection 
of the spacetime Einstein equations. Next, in Section 3, we discuss the relationship 
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between the constraints and evolution. Here, after first reviewing the proof of 
well-posedness, we discuss the ADM Hamiltonian formulation of Einstein's theory, 
noting the relationship between the Hamiltonian functional and the constraints. We 
also examine the evolution equations for the constraints, noting that if a set of data 
satisfies the constraints initially, it will continue to satisfy them for as long as the 
evolution continues. 

We discuss methods for constructing solutions of the Einstein constraint equa- 
tions in Section 4. We focus first on the most useful approach to date: the conformal 
method. After describing how the conformal method works, we discuss its success in 
parameterizing the set of all constant mean curvature (CMC) solutions of the con- 
straints, the difficulties which arise when constructing solutions with non-constant 
mean curvature, and the issue of finding physically relevant sets of data which satisfy 
the constraints. 

Closely related to the conformal method is the conformal thin sandwich ap- 
proach. We describe how it relates to the standard conformal method, and its 
major advantages and disadvantages. 

Instead of exploiting conformal variations to enforce the constraints, the original 
thin sandwich formulation |S1 12] varies the lapse (via an algebraic relation) and 
the shift vector. Under certain conditions this procedure leads to an elliptic system 
for the shift JJJ, an d an implicit function theorem then shows that the system 
is solvable for all nearby data. This constructs an open set of solutions of the 
constraint equations, from unconstrained data. However, the restrictions arising 
from the surprising cllipticity condition, that 7r y be positive or negative definite, 
mean that the original thin sandwich approach is viable only for a limited range 
of geometries, and at this stage it must be considered more of a curiosity than a 
practical solution technique. 

It is also possible to construct solutions of the Hamiltonian constraint by solving 
a semi-linear parabolic equation. Recent work 1711 l?2"] has shown that the es- 
sential feature of the original "quasi-spherical" construction ^1] can be generalized 
beyond the quasi-spherical foliation condition to give a flexible technique for con- 
structing exterior metrics of prescribed scalar curvature, satisfying geometric inner 
boundary conditions. Whilst the resulting parabolic method is restricted by the 
requirement that the background must be quasi-convex, it is able to handle several 
questions involving the constraints which cannot be addressed using the conformal 
method. 

We finish our review of constraint construction techniques by discussing the 
idea of gluing, and how this idea has been implemented to date. Two approaches 
to gluing have been developed and applied to solutions of the constraint equations. 
Connected sum or IMP gluing [22J builds new solutions by adding a cylindrical 
bridge (or wormhole) connecting a pair of points either on a single given solution or 
on a pair of given solutions. The Corvino-Schoen technique jSHlEZIESI gives a local 
projection from approximate solutions to exact solutions of the constraints. This 
may then be applied to smoothly attach a finite region in a given asymptotically 
Euclidean solution to an exterior Schwarzschild or Kerr solution, with a transition 
region connecting the two regions. Both approaches have been very useful in an- 
swering a number of longstanding questions regarding solutions of the constraints. 

We conclude this paper with a number of comments on important issues con- 
cerning the constraints which need to be addressed. 

2 Deriving the Constraints 

The Cauchy formulation is used primarily to construct new solutions of the Einstein 
gravitational field equations from specified initial data. The best way to understand 
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the origins of the constraint equations is to assume that we have a spacetime solution 
and to consider the induced data on spacelike hypersurfaces. We do this here, 
showing that the constraint equations necessarily must be satisfied by initial data 
(M 3 ,7, K), where M 3 is a manifold, 7 is a Riemannian metric on M 3 and if is a 
symmetric tensor on M 3 , if this data is to be induced by a spacelike hypersurfacc 
in a spacetime solution of Einstein's equations. 

2.1 3 + 1 Foliations of Spacetimes 

Let V 4 be a smooth 4-dimcnsional Lorentzian spacetime, with the smooth 1 metric 
g having signature (—1, 1, 1, 1). A hypersurface in V is an (embedded) submanifold 
M 3 V 4 of codimcnsion 1. M is spacelike if the induced bilinear form 7 := i*g 
is a Riemannian metric on M. Equivalently, M is spacelikc if at each point x G M 
there is a timclikc future unit normal vector n. We have the familiar diagram: 




If X, Y are any vector fields tangent to M, then as a consequence of the embed- 
ding, we can consider them as vectors in V and we can decompose the F-covariant 
derivative DxY into tangential and normal components, 

D x Y = V x Y + K(X,Y)n, (1) 

where V is the induced connection on M (it is also the Levi-Civita connection of 
the induced Riemannian metric on M), and if is a bilinear form (rank 2 tensor) 
on M, called the second fundamental form, or extrinsic curvature in the physics 
literature. From the relation 

K(X,Y)=g(D x n,Y) (2) 

and from the fact that the Lie bracket of X and Y is tangent to M, [X, Y] € TM, 
we find that if is a symmetric form, K (X, Y) = K(Y, X). 

A function t € C 1 (F) is a time function if its gradient is everywhere timelike. 
We say that a time function t is adapted to the hypersurface M if M is a level set 
of t, in which case we can choose adapted local coordinates (x, t). In terms of such 
coordinates the normal vector field takes the familiar lapse-shift form 

n = N- 1 (d t -X l d l ) , (3) 

where N is the lapse function and X = X t di is the shift vector. Equivalently, 

d t = Nn + X. (4) 

This time evolution vector field dt need not necessarily be timelike everywhere. (The 
shift vector X e TM is of course necessarily spacelike wherever it is nonvanishing.) 
However, if dt is timelike then the x = const, paths together make up a timclikc 
congruence of spacetime observers. 

Although it is useful in addressing certain questions to consider lower regularity, for our present 
purposes it is most convenient to assume that the spacetime and the metric are C°° . 
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The spacetime metric can be expressed in terms of the lapse and the shift and the 
spatial metric in adapted coordinates by the formula 2 

g = -N 2 dt 2 + (dx l + X 1 dt)(dx j + X j dt); (5) 

and the second fundamental form is given by 

K {j = K(d h 8,) = \N-\d tlij - C xllJ ) (6) 

where 

Cxla = Vi-Xj- + VjXi (7) 

is the Lie derivative in M of the spatial metric 7 in the direction X £ TM. Thus 
we obtain an expression for the d t evolution of the spatial metric, 

da ij =2NK ij +£ x j ij . (8) 

2.2 The Gauss and Codazzi Equations and the Constraints 

Based on the orthogonal decomposition of the covariant derivative operator D = 
V + K, and on the definition of the curvature on V, the Gauss and Codazzi equa- 
tions relate certain spacetime curvature components to spatial curvature and other 
expressions formed solely from data intrinsic to the submanifold M. For example, 
starting from the expression for the curvature on V 

Riem v (X, Y, Z, W) = g((D x D Y - D Y D X - D [X . Y] )Z, W), (9) 

if X, Y, Z, W are all vector fields tangent to M, then we find by a simple computation 
that 

Riem v (X,Y,Z, W) 

= Riem M (X, Y, Z, W) + K(X, W)K(Y, Z) - K(X, Z)K(Y,W). (10) 

This is the Gauss equation; it shows that the intrinsic curvature of M, measured 
by Riem M , is determined by curvature components of the ambient spacetime V and 
the second fundamental form K. In the perhaps more familiar index notation, the 
Gauss equation takes the form 

R&i = R&i + K a K jk - K ik K jU (11) 

where the indices k, I refer to an (unspecified) basis of the spatial tangent space 
TM, which may be equally well determined either by a coordinate basis di (holo- 
nomic basis) or by an orthonormal basis ej, i = 1, . . . , 3. 
A similar computation gives the Codazzi identity, 

Riem v (X,Y,n,Z) = g{(D x D Y - D Y D X - D [XiY] )n, Z) 
= V X K(Y,Z)- W Y K(X,Z) , (12) 

2 Our index conventions here are: mid latin alphabet indices etc.) run from 1 to 3 and 

correspond to directions tangent to the spacelike manifold M, while early latin indices (a, b, etc.) 
run from to 3 and correspond to spacetime directions. 
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where 

V X K(Y,Z) = D X (K(Y,Z))-K(V X Y 1 Z)-K(Y 1 V X Z) (13) 

is the covariant derivative in M of the tensor K. In index notation the Codazzi 
identity takes the form 

What about Riem v (n, Y, n, Z)l A calculation similar to those used to derive the 
Gauss and Codazzi equations does not lead to an expression solely in terms of 7, 
K, and their spatial derivatives. If we introduce the Lie derivative 

C n K(Y, Z) = D n (K(Y, Z)) - K(C n Y, Z) - K(Y, C n Z), 

along with the lapse and its spatial derivatives, then we have 

Riem v (Y, n, n, Z) = -£ n K(Y, Z) + K 2 (Y, Z) + iV- 1 Hess(7V)(y, Z), (14) 

where Kfj := K^Kjk and Hess(A r ) = V 2 iV is the Hessian or second covariant 
derivative of the lapse function N. This is sometimes referred to as the Mainardi 
equation, and has the index form 

R*mi = -UK i3 + K*K jk + N^Nm . (15) 

The index expressions for the Gauss, Codazzi and Mainardi equations given above 
are useful for calculating the 3 + 1 expressions for contractions of the Riemann 
tensor, such as the Ricci curvature tensor, the scalar curvature, and the Einstein 
tensor. For example, the Einstein tensor G a b — Ric^ — ^R y g a b satisfies 

2G„„ = R M + (tr 1 K) 2 - \\K\\ 2 (16) 
Gm = Klj-K^. (17) 

We note that these expressions for G nn and Gi n involve just 7, K, and their spatial 
derivatives. Hence it follows from the Einstein field equations 

Gab — 87rT ah (18) 

that the G nn and Gi n equations are constraints on the choice of the initial data 7 
and K . These are the Einstein constraint equations 

R M + (tr 7 X) 2 -||X|| 2 = 167rT n „, (19) 
Klj-K^ = 8nT ni . (20) 

Here T ao is the stress-energy tensor and describes the matter content of the ambient 
spacetime. Although T a b will itself generally satisfy some evolution equation, we 
will regard it as a prescribed field and often focus for definiteness on the case T a b = 
of vacuum (matter-free) spacetimes. 

It is clear from our discussion here that the constraint equations are conditions 
on the data (Af, 7, K) which are necessary for the data to arise from a spacelike 
hypersurface in a spacetime satisfying the Einstein field equations Ijl8|l . That these 
conditions are also sufficient (for the existence of a spacetime satisfying the Einstein 
equations and which induces the data on some hypersurface) is the content of the 
fundamental theorem of Y. Choquct-Bruhat 43 , which we discuss in the next 
section. 
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3 The Constraints and Evolution 



The constraint equations comprise four of the ten Einstein equations. The remaining 
six equations 

G jk = totT jk (21) 

describe how the data (7, K) evolve in a spacetime solution. This can be seen by 
rewriting the curvature terms in l|14|) . giving 

C„h ls = Gij - \G\ u,j + 2K?j - tr. hh, s - liic,;/ + N^ 1 N-ij . (22) 

Yet, as noted in the introduction, the constraints play a role in the evolution as 
well. We discuss this role here. We first show that the constraints are sufficient as 
well as necessary for a spacetime solution to evolve from a given set of data. We 
next explore the role of the constraints in the canonical, Hamiltonian, formulation 
of Einstein's theory. Finally, we discuss the evolution of the constraint functions 
under the evolution of the data generated by i|21|) . 

3.1 Well-posedness of the Vacuum Einstein Equations 

The Cauchy formulation of a given field theory is well-posed if (a) for any choice of 
initial data of specified regularity, there exists a solution which is consistent with 
that data, and (b) the map from the space of initial data to solutions is continuous. 
As shown by Choquet-Bruhat 03| 1 the Cauchy formulation of Einstein's (vacuum) 
theory of gravitation, for smooth initial data, is well-posed in coordinates of wave 
map gauge (harmonic) type. We now review the proof of this result 1431 EHl 1401 

Hid]. 

The vacuum Einstein equations take the form Kic v (g) = 0. If this were a 
hyperbolic system, well-posedness would follow from general results on quasi-linear 
hyperbolic systems (eg. jHHl El)- However, the system is not hyperbolic, which 
we verify by noting that by diffcomorphism invariance we have $*(Ric y ( 3 )) = 
Ric v ($*(ff)) for any diffcomorphism $ of V, so the symbol (leading order terms in 
the linearisation) of Kic v (g), considered as a partial differential equation in local 
coordinates, has a very large kernel. 

To get around this difficulty, we consider the reduced Einstein equations I39| 

p(g) := mc v (g) + \H w g = 0, (23) 

where Cw9ab = D a Wb + D]yW a is the Lie derivative of the metric g in the direction 
of the vector field W, and W is chosen to have leading terms 

W a = -g cd (d c g da - \d a g cd ) + F{g) (24) 

where the terms F(g) do not involve derivatives of g. A short calculation shows 
that p{g) has symbol ^g oh £ a £b, exactly that of the wave equation of the Lorentzian 
metric g itself. Thus p(g) = forms a quasi-linear hyperbolic system, 

9 cd d 2 cd g ab = Qa b (g,dg), (25) 

where Q a b is quadratic in dg and depends also on the F(g) terms in (I24|) . As noted 
above, the initial value problem for p{g) = is well-posed. 

However in general solutions of p{g) = will not satisfy the vacuum Einstein 
equations. The idea now is to choose the initial conditions (and, perhaps, the 
precise form of the terms F(g)) in such a way that the vector field W can be shown 
to vanish identically 
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To understand the meaning of W , recall the definition of the tension field u)((j>) 
of a map 0: (V,g)^(N,V), 

w a (0 = V a D a( j) a = g cd (d 2 cd r - T e cd d e r + r^dcfd^) , (26) 

where N is some chosen four dimensional manifold and V is an arbitrary connection 
on N. If we now assume is a diffeomorphism, which we use to identify the 
coordinates x a on V with y a on N by pullback x a — <fi*(y a ), then the tension 
becomes 

w <»(0) = -g ab g cd (V c9db - ±V bffcd ) , (27) 

which shows by comparison with 124J1 that W may be chosen as the metric dual of 
the push- forward by </> _1 of the tension, 

W a = g ab <f>Z 1 {u b {4>)) . (28) 

Note that the reduction to a quasi-linear hyperbolic system holds for any choice 
of the target (N, V), where the connection V need not be compatible with any 
metric on N. In the simplest case V is the flat connection on IR 4 and the condition 
W = translates into O g x a = in local coordinates x a on V and N. The gauge 
condition W(<f>) — is known as the wave map or Lorentz-harmonic coordinate 
gauge. The common terminology "harmonic gauge" is inappropriate since solutions 
of the wave equation do not have the mean- value ( "harmonic" ) property implied by 
the corresponding Euclidean Laplace equation. 

Good existence theorems for the Cauchy initial value problem for the reduced 
Einstein system (|25|) are well-known (see, for example, [741 1581 16*5] ). and the initial 
data (3ah(0), d t g a b(0)) are freely prescribable, subject only to the condition that 
5ah(0) has Lorentz signature with Mq spacelike. In particular, if the geometric data 
(jij,Kij) are given (subject only to the vacuum constraint equations i|19 fl -l2U |l 1. 
then for any choice of the lapse and shift g$ a we may recover dtgij via (JSJ. Since 
9ij = Jijj w e see that initial data (g o h(0), dtg a b{0)) may be chosen consistent with 
(7, if) for any choice of (fifoo(O), d t g 0a (0)). 

Of course, it is not sufficient to just solve the reduced Einstein equations; in 
order to recover a solution of the full vacuum Einstein equations we must also 
ensure that the solution satisfies W = 0. However, it is an easy consequence of the 
second Bianchi identities and p(g) = that W must satisfy 

□ g W + Ric v \g)W = 0. (29) 

Thus we aim to construct the full initial data (<7ab(0), dtgab(0)) for the evolution 
equation p(g) = in such a way that W(0) = and dtW(0) = 0. If we can do this, 
then uniqueness for solutions of the Cauchy problem for the hyperbolic evolution 
equation Q29J1 will imply that W = everywhere. 

To ensure that VF(0) = we may either choose the background metric h on 
N appropriately or we may expand Q24J) and choose the components dtgoa(0) ap- 
propriately. As noted above, the choice of dtgo a (0) does not affect the prescribed 
geometric data (7, if). 

To find the conditions necessary for choosing initial data with ^^(O) we com- 
bine the reduced Einstein equation I123H and the condition W(0) = to obtain 
— 2Ric^ fc = d a Wb + dbW a at t = 0, which may be rearranged into the form 

g 00 d t W a = -2G° + (terms vanishing at t = 0). (30) 

Then, noting that G° a are the constraint functions, we see from l|30() that the con- 
dition that the geometric data (7, if) satisfy the constraint equations ensures that 
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dtW(0) — 0. Since solutions to l|29(l are uniquely determined by the initial data 
(W(0), d t W(0)), the corresponding solution of the reduced equation iffify will have 
W = everywhere and thus g will satisfy the vacuum Einstein equations with geo- 
metric data (7, A"). This shows that the constraint equations are also sufficient for 
the existence of a compatible vacuum spacetime, as claimed. 

In this formulation of the Einstein equations, the lapse and shift are determined 
from goa(t) by solving the reduced Einstein equations. Alternatively, we may ask 
whether it is possible to prescribe the lapse and shift as given functions of x and 
t (for example, N(x,t) — 1 and X(x,t) = for all x <E M and t > 0). As shown 
recently by Choquet-Bruhat [28] , this is possible although in a weaker sense than the 
Einstein system in the wave map gauge described above. Specifically, the Einstein 
system for general lapse and shift is "non-strictly hyperbolic"; as a consequence, 
the system is well-posed in Gevrey spaces, and so has local existence. We note 
as well that a number of researchers [311 EH] have formulated extended first order 
versions of Einstein's vacuum equations and have verified symmetric hyperbolicity 
and therefore well-posedness in more standard spaces for these versions. 

Because the initial lapse-shift variables goa(0) (cf. 10) are freely specifiable in the 
above construction, there remains the question of the uniqueness of the spacetime 
development of a given set of initial data. It turns out that solutions constructed 
using different choices of the lapse and the shift but the same geometric data (7, K) 
are in fact equal in some neighbourhood of the initial surface, after suitable coordi- 
nate changes |29) . This property is established by constructing a wave map ip from 
one solution to the other, and using the identification of the initial data to ensure 
the wave map initially at least is a spacetime isometry. The two metrics can then 
be compared directly (via g\ and ip*(g2)), and uniqueness for the reduced Einstein 
equations leads to geometric uniqueness. This shows that the geometric data (7, K) 
determine the resulting spacetime uniquely, at least locally. 

3.2 Einstein's Theory as a Dynamical System and the ADM 
Hamiltonian Formulation 

If we fix a choice of the lapse and shift for which the system of Einstein's equations 
is well-posed (see the discussion above), then Einstein's equations can be viewed for- 
mally as a dynamical system. More specifically, fixing a three dimensional manifold 
M and letting TM. denote the tangent space to the space M of smooth Riemannian 
metrics on M, we find that (with a fixed choice of lapse N and shift X) the Einstein 
evolution equations 

2NK ij + C x j i:i , (31) 

VijN + CxK i:j + N(RicYj - Ricf/ + 2K?K jk - tr 7 X K l3 ) (32) 

specify a flow on TM . The Kij evolution equation follows from the Mainardi equa- 
tion, cf. 114|l and l[22[l. For the vacuum Einstein case Ric^ = 0; in the nonvacuum 
case, Ric]^ depends on the matter fields through the stress-energy tensor T a \>. 

There is an important sense in which the constraints generate the flow, which 
we see by recasting Einstein's theory in the canonical Hamiltonian form. Originally 
motivated by attempts to obtain a theory of gravity consistent with quantum theory, 
and often labeled the "ADM" Hamiltonian formulation because of the early work 
on these ideas by Arnowitt, Deser and Misner [7], this is based on the analysis of 
the Einstein-Hilbert Lagrangian J v R v dv g . 

One way to approach this analysis is via the second variation formula for space- 
like hypersurfaces [701 IIP): for more traditional treatments see [SI ESI EZ]- We 



dt 
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let F : (— e, e) X M — » y be a 1-parameter family of embeddings s : M i— > M s = 
F(s,M), with variation vector y = -F*(Jj), with normal vector n — n s and with 
second fundamental form K = K s along M s . If y is everywhere transverse to M 
then s i— > M s is a foliation and we can regard Y, n, and K as fields on V; otherwise 
they live naturally on (— e, e) x M via pullback. The first variation of the area \M S \ 
is given by 

Dy{dv 1 ) = —g{Y,n)tv 1 K dv-y , (33) 
and the second variation formula is 



Dyt^K = g(Y, n)(\\K\\ 2 + Ric y (n, n)) + g(Y, Vtr 7 K) - A M g(Y, n) , (34) 

where Dy = and the s-dependence of n and K is understood. This follows also 
from the Mainardi equation 1)150. Since from IjlUI) and 1(12(1 we have 



2Ric y (n, n) 



R v 



R M - \\K\\' Z + (t^Kf 



and since — g(Y, n) — N is the lapse if y = is the evolution vector of a foliation 
by spacelike hypersurfaces, we find that 134(1 can be rewritten as 



R 



R 



|K|| 2 + (tr 7 if) 2 - 2N~ 1 A M N + 2D n (tTjK) 



(35) 



Integrating this quantity over the region bounded by two spacelike hypersurfaces 
F({sq} x M) and F({si} x M), we obtain the Lagrangian 



R v dv v 

F([«o,«i]xM) 
si 

R M + \\K\\ 2 - (tr y Kf) N - 2A M Nd Vl ds 



2tr 7 if dv^ , 



(36) 



using 1(33(1 . From this Lagrangian, we calculate ir^ , the momentum conjugate to 
7 y (see |Z5]) as 

<5£ 



= v/det7 (iT J - tr 7 K7 l 



(37) 



Then ignoring boundary terms, we calculate the Hamiltonian TL = J M ^ lJA /ij — £ in 
ADM form: 



Ha 



DM 



(38) 



M 



where £ a = (N, X 1 ) is the four- vector consisting of the lapse and the shift, and 
where 



$o = R M det(7) 1/2 -dct(7)- 1 / 2 (||^|| 2 -i(tr 7 ^) 2 ) 
*< = - 27yV fc 7r jfe 



(39) 
(40) 



are the constraint operators, written in densitized form, in terms of the canonical 
variables (7, it) rather than (7, K). It is easily verified that the equations $ a (7, K) — 
0, a = 0, . . . , 3 are equivalent to (HHJ,© with T 0a = 0. 

Given an explicit expression for the Hamiltonian such as (|38H . we can readily 
calculate the evolution equations for (7, ir) and for any functional of these quantities. 
In particular, we have the ADM form [J] of the evolution equations 1)31(1 . 1(32(1 



d 


7 




1 " 


It 


7T 




-1 



D* ( 7> 7r)*£°, 



(41) 
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where DQ* is the formal adjoint of the linearization (or functional derivative oper- 
ator) £>$ a . Explicitly we have [4^] : 

D$(7,7r)(/i,p) = 

((SSh - Atr 7 /i - fty (Ric My - ±R M 7 y ')) ^7 

+ (tr 7 7T7T« - 2tt^ + i|7r| 2 7 « - £(tr 7 7r)V 3 ') /^7 
+ p lJ (tr 7 7T7ij - 27Ty) /\/7 , 

7r jfe (2V J /i 2fc - Vjhj-fc) + 2^V fe7 r jfe + 2 7lfe V,p J ^ , (42) 

D$(7,7r)*C - 

f (V*V J JV - AN-f ij - /VRic Mij + |iVR M 7 ij ') ^7 

+ TV (tr 7 7T7r y ' - 24^^' + i|7r| 2 7 y ' - ±(tr 7 7r)V) / yft 

+ (x k V k n ij + V k X k n ij - 2V fc A (l 7! J ' ) ' £ ) , 

N (tr 7 7r 7ij - 2tt«) /V7 - 2^V (i X i) ) , (43) 

where SSh = V l V J % and £ = (N,X l ). 

Note that in the literature there are a number of different explicit expressions 
for the time derivatives of 7 and it [3]. These expressions can all be related by the 
addition of terms which vanish if the constraint equations are satisfied. 

The adjoint used in l|41f) is formal, in that it ignores the contributions of bound- 
ary terms. For spatially compact spacetimes, there are of course no boundary terms. 
For asymptotically flat spacetimes, the resulting asymptotic boundary terms in the 
Hamiltonian are closely related to the ADM total energy-momentum jJSI , and they 
play a major role in the analysis of the physics of such spacetimes. 

An important consequence of the form 1)41(1 of the evolution equations is the re- 
sult of Moncrief |63| . which shows that a vector field £ = (N, X 1 ) satisfies -D$( 7 , K)*£ = 
for sufficiently smooth vacuum initial data (7, K) if and only if £ is the restriction 
to M of a Killing field in the spacetime development of the initial data. 

3.3 Preserving the Constraints 

If (70,710) is initial data which satisfies the (vacuum) constraint equations and has 
vacuum evolution (-f(t),w(t)) corresponding to some choice of lapse-shift (N, X 1 ), 
it is natural to ask if the evolving data satisfies the constraints for all values of 
t. A first step in showing this is to calculate the time derivative of the constraint 
operators H = $ and Ji = cf. 0: 

X7 m (X m H) + 2-yVViNJj + NV l J i + \N 1 - 1 '^H lij % ij , (44) 
2WiNH + NW l H + Vt{X k Ji) + W l X k J k . (45) 

These formulas follow directly from the conservation law V b G a & = and the fact 
that the evolution equation l(32|l implies GV, =0, 1 < i,j < 3. The virtue of these 
expressions for the time derivatives of the constraints is that together they comprise 
a symmetric hyperbolic system. We may then use standard energy arguments (based 
on an energy quantity which is quadratic in H and Ji) to argue that if initially H and 
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Ij. 
dt l 
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Ji vanish, then for all time consistent with a foliation of the spacetime development 
of the initial data, H and Jj vanish as well. 

The formal Hamiltonian interpretation of the evolution equations (|41J) regards 
(7(f), 7r(t)) as defining a trajectory in the cotangent bundle T* M. of the space M. 
of Riemannian metrics 7^ on M. It is straightforward to endow this phase space 
T with the structure of a Hilbert or Banach manifold, for example ('jij,^ 3 ) € 
H S (M) x H S ~ 1 (M), s > 2 for compact M using the Sobolev spaces H s , and then 
the evolution (|41|l is determined from a densely defined vector field on T . The fact 
that the constraints are preserved shows that a trajectory of the evolution (|41(1 lies 
in the constraint set C, the subset of J- satisfying the constraint equations, provided 
the initial data (7(0),tt(0)) lies in C. 

Note that while the constraints are preserved under the explicit, smooth, evolu- 
tion generated by (|41|l . if we evolve instead numerically (with the consequent un- 
avoidable numerical errors), then the dynamical trajectory will leave the constraint 
set C. Indeed, numerical tests suggest that such a trajectory will exponentially 
diverge from C. This property, and how to control it, is under intensive study by 
numerical relativists; for a recent example see 

3.4 Linearization Stability 

The linearized Einstein equations govern the evolution of perturbations 5g of a 
given background spacetime metric go and have been extensively analysed, in view 
of their importance in astrophysics and cosmology. It is then natural to ask whether 
linearized solutions actually correspond to solutions to the full nonlinear equations; 
that is, whether a given solution of the linearized equations is tangent to some 
curve of solutions of the full equations. Such a perturbation solution is called 
integrable. By the fundamental local existence theorem |4H| of Yvonne Choquet- 
Bruhat, a perturbation is integrable if all solutions of the linearized constraint 
equations _D$( 7-7r )(/i,p) = correspond to full solutions of the constraints; i.e., if 
the constraint set C is a Hilbert submanifold. This is the question of linearization 
stability and it is somewhat surprising that this fails in certain situations. 

The implicit function theorem (see e.g. PP) shows that C is a Hilbert submanifold 
of J- provided the linearization Z)$( 7(7r j has closed range and is surjective. Now the 
range of £)$( 7iJr ) is i 2 (M)-orthogonal to the kernel of the adjoint -D$* 7 T y which 
corresponds by work of Moncrief E2J to space-time Killing vectors. Using this 
insight and elliptic theory for the operator D$( 7>7r )D$? it has been shown by 
Fischer and Marsden 021 EEH IHj that if the data (7, ir) does not admit Killing vectors 
and if M is compact then C is locally a Hilbert manifold, at least near smoother 
data. 

Moreover, at points of C corresponding to data for spacetimes with a Killing 
vector field, C is not a submanifold but instead instead has a cone-like singularity 
arising from a quadratic relation |41j . This condition arises from the second deriva- 
tive of a curve in C, A 1— > $(7(A), 7r(A), about (70,710) admitting a Killing field £. 
Integrating over M, discarding boundary terms and using D<$>f = gives 



which is the required quadratic condition on the constrained variations (h,p). In |64| 
it is shown that this condition is equivalent to the vanishing of the Taub quantity, 
constructed from a solution of the linearized Einstein equations about a spacetime 
admitting a Killing vector. 

The situation for asymptotically flat spacetimes differs significantly from the 
case of compact M. Using weighted Sobolev spaces (with local regularity 7 € 




(46) 
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H s , s = 2) it is possible to show that the C is everywhere a Hilbert submanifold 
although Killing vectors such as rotations are possible, they do not satisfy 
the asymptotic conditions needed to prevent surjectivity of D<I>( 7i7r ) and thus the 
quadratic conditions on constrained variations do not arise. This work [§] also 
shows that the asymptotically flat phase space T = (g + -ff£ 1 / 2 ) x P rov ides a 

natural setting in which the ADM Hamiltonian and energy-momentum functionals 
become smooth functions, for example. 

4 Solving the Constraint Equations 

There are three reasons for seeking to construct and study solutions of the Einstein 
constraint equations. First, we would like to obtain initial data sets which model the 
initial states of physical systems. Evolving such data, we can model the gravitational 
physics of those physical systems. Second, we would like to understand the space 
of all globally hyperbolic solutions of Einstein's gravitational field equations. Since, 
as noted earlier, Einstein's equations are well-posed, the space of solutions of the 
constraints parametrizes the space of solutions of the spacetime field equations. 
Third, we would like to know enough about the nature of the space of solutions of 
the constraints and the Hamiltonian dynamics of the classical Einstein equations 
on this space to be able to consider descriptions of gravitational physics which are 
consistent with the quantum principle (i.e., to "quantize gravity"). 

The easiest way to find a set of initial data (M 3 , 7, K) which satisfies the Einstein 
constraint equations is to first find an explicit solution (V , g) of the spacetime 
Einstein equations (e.g., the Schwarzschild or Kerr solutions), and then to choose 
a Cauchy surface M 3 °-> V 4 in that spacetime solution: The induced Riemannian 
metric 7 and second fundamental form K on M 3 together solve the constraints. 
Since the set of known spacetime solutions is very limited, this procedure is not 
especially useful for making progress towards the goals discussed above (although 
we shall see below that solutions of the constraints obtained from slices of known 
spacetime solutions can be very handy as building blocks for "gluing constructions" 
of solutions of the constraints, which are in turn potentially very useful for modeling 
physical systems). 

Our focus here is on more comprehensive methods for finding and studying 
solutions of the constraints. We first discuss the conformal method. This is by far 
the most widely used approach, both for physical modeling and for mathematical 
analyses, and so it receives the bulk of our attention here. We next consider the 
thin sandwich and conformal thin sandwich ideas. The latter is closely related 
to the conformal method; we compare and contrast the two here. The fourth 
approach we consider is the quasi-spherical ansatz. Although not as comprehensive 
as the conformal method, the quasi-spherical ansatz is potentially useful for certain 
key applications, which we discuss below. We conclude by describing the recently 
developed procedures for gluing together known solutions of the constraints, thereby 
producing new ones. 

For convenience, we work primarily in this paper with the vacuum constraint 
equations. We note, however, that most of what we discuss here generalizes to the 
Einstein-Maxwell, Einstein- Yang-Mills, Einstein-fluid, and other such nonvacuum 
field theories. 

4.1 The Conformal Method 

The vacuum Einstein constraint equations (|19|l . (|20|l with T nn = and T n i = 
constitute an underdetermined system of four equations to be solved for twelve 
unknowns 7 and K. The idea of the conformal method is to divide the initial data 
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on M 3 into two sets — the "Free (Conformal) Data", and the "Determined Data" — 
in such a way that, given a choice of the free data, the constraint equations become 
a determined elliptic PDE system to be solved for the Determined Data. 

There are at least two ways to do this (see [35] for others). The sets of Free 
Data and Determined Data are the same for both procedures; we have 

Free ("Conformal") Data: 

Xij a Riemannian metric, specified up to conformal factor; 

<Tij a divergence- free 3 (Wy = 0), trace-free (A y <7jj = 0) 

symmetric tensor; 
r — a scalar field, 

Determined Data: 

4> a positive definite scalar field, 

W % — a vector field. 

The difference between the two procedures has to do with the form of the equations 
to be solved for the Determined Data, and the way in which the two sets of data are 
combined to obtain 7 and K satisfying the constraints. In one of the procedures, 
which we label the "Semi-Decoupling Split" (historically called "Method A"), the 
equations for W and <j> take the form 

Semi-Decoupling Split 

V^LW)) = |0 6 V,r (47) 
A0 = \R<t>-\{a ij +LW ij ){a ij +LW ij )<p- 7 + ±i*<l> h , (48) 

where the Laplacian A and the scalar curvature R are based on the Ay -compatible 
covariant derivative Vj, where L is the corresponding conformal Killing operator, 
defined by 

(LW)ij = ViWj + VjWi - fAyVfcW^, (49) 
and we construct 7 and K from the free and the determined data as follows 

lij = 0% (50) 
Ka = ^(aij+LWi^ + ^XijT. (51) 

Using the other procedure, which we label the "Conformally Covariant Split" (his- 
torically "Method B"), the equations for W and <fi are 



Conformally Covariant Split 



V l {LW)) = -§{LW))VM<i> (52) 

+ ±(r 2 ~(LWr(LW) lJ )cb 5 (53) 

and the formulas for 7 and K are 

lij = 4 A. (J (54) 
K tJ = ^aij+^LWij + ^XijT. (55) 

Each of these two methods has certain advantages. For the semi-decoupling split, 
if one chooses the mean curvature r to be constant, then the <$> dependence drops 

3 In the free data, the divergence-free condition is defined using the Levi-Civita covariant deriva- 
tive compatible with the conformal metric Xij. 
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from (|47J) . and the the focus of the analysis is on (|48J) . the Lichnerowicz equation. 
This is not true for the conformally covariant split; however in this latter case, one 
finds that a solution exists for free data (Ay, cn 3 r) if and only if a solution exists 
for (8 4 Xij, 8~ 2 (iij, t). Far more is known mathematically about the semi-decoupling 
split, and this approach has been used much more in applications, so we focus on 
the semi-decoupling split in the rest of this paper. We do note, however, that 
numerical relativists have very recently begun to apply the conformally covariant 
split for certain studies jHSJ- 

Is it true that, for every choice of the free data (Ay , <7ij , t) , one can always 
solve equations l|17|) and (|4*5|) for the determined data (<p, W) and thereby obtain 
a solution of the constraint equations? It is easy to see that this is not the case: 
Let us choose the manifold M 3 to be the three sphere, and on M 3 we choose a 
metric A with non-negative scalar curvature, we choose a to be zero everywhere, 
and we choose r to be unity everywhere. One readily verifies that every solution 
to the equation \/i(LW)j — has LWy = 0. The Lichnerowicz equation then 
becomes = R<f) + <j) 5 . Since the right hand side of this equation is positive 
definite (recall the requirement that > 0), it follows from the maximum principle 
on closed (compact without boundary) manifolds that there is no solution. 

In light of this example, the main question is really the following: For which 
choices of the manifold M 3 and the free data can one solve l|47|) and (|48|l ? It turns 
out that we know quite a bit about the answer to this question, yet still have quite 
a bit to learn as well. To describe what we know and do not know, it is useful to 
categorize the question using the following criteria: 

Manifold and Asymptotic Conditions 

• Closed 

• Asymptotically Euclidean 

• Asymptotically hyperbolic 

• Compact with boundary conditions 

• Asymptotically Euclidean with interior boundary conditions 
Regularity 

• Analytic 

• Smooth 

• Sobolev and Holder classes 

• Weak 

• Asymptotic fall off conditions at infinity 
Metric Conformal Classes 

• Yamabe positive 

• Yamabe zero 

• Yamabe negative 

Mean Curvature 

• Constant ("CMC") 

• Near Constant (small max l^ll ) 

^ min |t| ' 

• Non Constant 
a Classes 
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• er is zero everywhere on M 3 

• a not zero everywhere on M 3 

The mean curvature of the data turns out to be the most important factor in 
separating those sets of free data for which we know whether or not a solution exists 
from those sets for which we do not. In fact, if the mean curvature is either constant 
or near constant, we can almost completely determine whether or not a solution 
exists (at least in those cases with no boundaries present). On the other hand, if 
the mean curvature is neither constant nor near constant, we know very little. For 
sufficiently smooth free data, we may summarize the known results as follows: 

Constant Mean Curvature 

• Closed: Completely determined. For all but certain special cases (see 
example discussed above), solutions exist. See [2] and |BT1] . 

• Asymptotically Euclidean: Completely determined. For all but certain 
special cases (depending entirely on the conformal class), solutions exist. 
See and El]. 

• Asymptotically Hyperbolic: Completely determined. Solutions always 
exist. See jH] and @] 

• Compact with Boundary: Some cases determined, most unresolved. 

• Asymptotically Euclidean with Interior Boundary: Some cases deter- 
mined, some unresolved. See |BJ and 

Near Constant Mean Curvature 

• Closed: Mostly determined. For all but certain cases, solutions exist. 
One special case known of nonexistence. A small number of special cases 
unresolved. See [HI], [S3 an d [EE] 

• Asymptotically Euclidean: Mostly determined. No cases known of nonex- 
istence. Some cases unresolved. See [30] 

• Asymptotically Hyperbolic: Mostly determined. No cases known of nonex- 
istence. A small number of cases unresolved. See |56| 

• Compact with Boundary: Small number of cases determined, most unre- 
solved. 

• Asymptotically Euclidean with Interior Boundary: Nothing determined. 
Non Constant Mean Curvature Very little known. 

This summary is the compilation of results from a large number of works (note 
references listed above), including some not yet written up. The summary is fairly 
sketchy, since there is inadequate space here to include many of the details. We 
do, however, wish to give a flavor of what these results say more precisely, and how 
they are proven. To do this we shall discuss a few representative sub-cases. 

We first consider the case in which the manifold is presumed closed, and in which 
we choose free data with constant r. As noted above, as a consequence of the CMC 
condition, the equation (|47(l for W is readily solved, and we have LW = 0. What 
remains is the Lichnerowicz equation, which takes the form 

= \R<j) - \a ij airf-l + j^T 2 (j) b . (56) 

There are three key PDE analysis theorems which allow us to determine, for 
any set of free data consistent with these conditions, whether or not a solution to 
the Lichnerowicz equation exists. We state each of these theorems in a form most 
suited for this purpose (See for much more detail): 
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1. The Maximum Principle: If we choose free data such that for any <f> > the 
right hand side of the Lichnerowicz equation Ij56|l is cither positive definite 
or negative definite, then (with M 3 presumed closed) for that free data there 
can be no solution. 

2. The Yamabe Theorem: For any given Riemannian metric A on a closed three 
manifold, there is a conformally related metric 4 X which has constant scalar 
curvature. For fixed A, all such conformally-related constant scalar curvature 
metrics have the same sign for the scalar curvature. This allows one to parti- 
tion the set of all Riemannian metrics on M 3 into three classes, depending on 
this sign. We label these Yamabe classes as follows: y + (M 3 ), y°(M 3 ), and 

y-(M 3 ). 

3. The Sub and Super Solution Theorem: If there exist a pair of positive functions 

> 4>— such that 

A0+ < ±i?</>+ - + ±t 2 4>1 (57) 

and 

A0_ > \R4>- - la ij a ij( f)Z 7 + t^t 2 ^, (58) 
then there exists a solution <f> of the Lichnerowicz equation (|56() , with <f> + > 

<t>><t>- 

The Yamabe Theorem is very useful because the Lichnerowicz equation in the 
form (|56Jl (with LW — 0) is conformally covariant in the sense that it has a solution 
for a set of free data (A, a, r) if and only if some function tp > it has a solution 
for the free data (V> A, ip~ 2 cr, r). Hence, to check solubility for a given set of free 
data, we can always first perform a preliminary conformal transformation on the 
data and therefore work with data for which the scalar curvature is constant. 

It is now straightforward to use the Maximum Principle (together with the 
Yamabe Theorem) to show that there is a collection of classes of CMC free data for 
which no solution exists. This holds for all of the following classes of free data: 

(A e y + (M 3 ), a = 0, t = 0), (A € y + (M 3 ), a = 0, r ^ 0), 

(A e y°(M 3 ), a = 0, r ^ 0), (A E y°(M 3 ), a 0, r = 0), 

(A e y- (A/ 3 ), a = 0, r = 0), {Xey- {M 3 ),a ? 0, r ^ 0). 

Note that here,"cr ^ 0" means that the tensor a is not identically zero on M 3 . 

It is true, though not so simple to show, that for all other classes of CMC free 
data on a closed manifold, a solution does exist. We show this using the sub and 
super solution theorem. For free data of the type (A S y~(M 3 ), a = 0, r ^ 0), it is 
relatively easy to show that there are constant sub and super solutions. We merely 
need to find a constant <fi + sufficiently large so that —<fr+ + j^t 2 ^ is everywhere 
positive, and a constant 0_ sufficiently small (yet positive) so that — + yjT 2 ^ 1 ! 
is everywhere negative. Constant sub and super solutions are also readily found for 
free data of the type (A e y-(M 3 ),a ^ 0, r ^ 0), or (A e y°(M 3 ), a = 0,r = 0). 
For the remaining types of free data, we need to find either a nonconstant sub 
solution or a nonconstant super solution. This takes a bit of work; we show how to 
do it in 150]. 

It is worth noting that the proof of the sub and super solution theorem (See |5()j ) 
is a constructive one: Starting with the super solution </> + , one solves a sequence 
of linear equations, and one then shows that the monotonic sequence of solutions 
4>+ ^ 4>i ^ 02 > • • • > <t>- converges to a solution of the Lichnerowicz equation. 
This constructibility can be useful for numerical relativity, as has been shown in 
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We next consider the case in which the free data is chosen to be asymptotically 
Euclidean (See, for example, for a definition of this property), with r = 0. 
For asymptotically Euclidean metrics, a Yamabe type theorem again plays a key 
role. Proven by Brill and Cantor [53], with a correction by Maxwell |6 1 1 . this result 
says that if A is asymptotically Euclidean, and if for every nonvanishing, compactly 
supported, smooth function / we have 



(,,(|V/| 2 + i?/ 2 )%/detA , 
inf ■ iA/VI ,. „ ; > 0, (59) 



{fm \\j\\ L 2, 

then for some conformal factor 0, the scalar curvature of 9 A is zero. Moreover, 
if an asymptotically Euclidean metric fails to satisfy this condition, then there is 
no such transformation. Metrics which do satisfy this condition have been labeled 
(somewhat misleadingly) as "positive Yamabe metrics" . 

The positive Yamabe property just defined is exactly the condition which deter- 
mines if, for a set of asymptotically Euclidean free data, the Lichnerowicz equation 
admits a solution. That is, as proven in |26| . (|56|) admits a solution <fi with suit- 
able asymptotic properties if and only if the (asymptotically Euclidean) free data 
(A, <j,t — 0) has positive Yamabe metric A. This is true regardless of whether a 
vanishes or not. This result is not proven directly using a sub and super solution 
theorem, but the proof does involve a converging sequence of solutions of linear 
equations, and is therefore again constructive. 

The last case we consider here is that of near constant mean curvature free data 
on closed manifolds. This case is more complicated than the CMC case, since now 
we must work with the coupled system (|47ll - H48|) . However, with sufficient control 
over the gradient of r, we can handle the coupled system, and determine whether 
or not solutions exist for almost all sets of near CMC free data. 

The key extra tool we use to carry out this analysis is the elliptic estimate for 
the "vector Laplacian" operator VjL( )*■ appearing on the left side of (|47|l . Such a 
result holds for any linear, elliptic, invertible operator [22] such as this one. 4 One 
has the Elliptic Estimate: ||W J '||fl- k+a < C||Vj(LW)'-||fl- k , where Hk is the Sobolev 
space of square integrable vector fields for which the first k derivatives are square 
integrable as well, where 1 1 | \jj k is the corresponding norm, and where C is a constant 
depending on the chosen geometry (M 3 ,Ay). Based on this estimate, together 
with the Sobolev embedding theorem |22| and standard integration inequalities, we 
obtain a pointwise estimate of the form 

\LW\ < Cmax/maxlVrl, (60) 

M 3 M 3 

where C is also a constant depending only on the geometry (M 3 , Ay). This pointwise 
inequality (|60|l is the crucial estimate which allows us to handle the coupled system, 
for sufficiently small |Vr|. 

To see how to prove that a solution exists in one of these near CMC cases, let us 
consider free data with A S y~ (M 3 ),r > 0, and with maXjv?3 7"^ sufficiently small. 5 
We claim that, for such free data, a solution exists. To show this, we consider the 
sequence of semi-decoupled PDE systems 

Vi{LW {n) )\ = f^-D^r (61) 

A<fr n) = lflfr n) -l (<^ + LWfyfa + LW Mij )^ + ^T 2 ^ n) .(62) 

4 The operator V;L( is invertible only if A has no conformal Killing field. However, the 
pointwise estimate I6UI we obtain for \LW\ can be derived even if A has a conformal Killing field, 
using a slightly more complicated elliptic estimate than the one presented here. 

5 The statement that " — . m3 -S (is) sufficiently small" may seem nonsensical, since this 

quantity is dimensional. However, as seen in the proof 1541 . the more precise statement of this 
condition involves the constant C, which has dimensions as well. 
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Choosing a suitable initializing value for <fio , we first show that a sequence {4>( n ) , Wr n ) ) 
of solutions to [j6H62Jl exists. For the decoupled equation ifBTJl . the existence of a 
solution Wi n ) follows from the invertibility of the vector Laplacian. For Q62[l. ex- 
istence follows from the sub and super solution theorem, since we readily find a 
sequence of constant sub and super solutions for this choice of free data. We next 
show that there are (positive) uniform upper and lower bounds for the sequence 
(0(„)). This is where we first need to use the estimate (|60l) for |LM^(„)|: Specifically, 
after arguing that there is a uniform (constant) upper bound for the </>(„)'s if we 
can find a positive constant £ such that 

C 3 > ^ jC 2 + ■ 3 a ^ + , 3 2 LWg } LW {n)ij , (63) 

we use the |LW(„)| estimates to show that it is sufficient to find a positive constant 
£ satisfying 



^3 > J z-2 , 3 _ M _ maxMsjVrp 3 
2 mm A /3 r z mm M 3 2 min M 3 



We immediately see that if, for fixed A and tr, we choose Cmaxjyfa | Vr| 2 /(2 minjv/3 t 2 ) 
sufficiently small, then such a £ exists. The existence of the uniform upper bound 
for <frr n \ follows. 

Establishing these uniform bounds is the crucial step in our proof that solutions 
exist for free data of the type we are discussing here, as well as for other classes 
of near CMC free data. Once we have these bounds, we can carry through a 
contraction mapping argument to show that the sequence {4>( n ) > W(n) ) converges 
(a bit more squeezing of |Vr| is needed to carry this out). We can then go on to 
use continuity arguments to show that the limit is a weak solution, and finally use 
bootstrap arguments to show that in fact the weak solution is a strong solution, of 
the desired smoothness. The details of this proof are presented in |54j . 

A very similar argument can be used to prove that equations l|47 |l -l|48 |l admit 
solutions for a number of classes of near CMC free data, including the following: 

1. A S y-(M 3 ),r < 0, and n ™>' 3 JJ T ] 1 sufficiently small; 

2. A G y + (M 3 ), a not identically zero, and ^*^ f3 J^ sufficiently small; 

3. A S y°(M 3 ), t nowhere zero, a not identically zero, and °™" 3 ^ Vr ^ ,_ 

sufficiently small. 

What about the other near CMC cases on closed M 3 ? Just recently, with Niall 
O'Murchadha, we have found for the first time a class of near CMC free data for 
which a solution does not exist. The data is of the following type: 

A S y+(M 3 ) U y°(M 3 ),a = 0, t nowhere 0, and max ^ 3 l V [l sumcie ntly small. 

mm M 3 \t\ 

The proof that no solution exists for free data of this type is a relatively simple 
application of the pointwise estimate i|tj(J|) for \LW\ discussed above, and the maxi- 
mum principle. The details are presented in |55j . As for near CMC data on a closed 
manifold which is not one of the types discussed here, nothing is yet known. Note 
that, in a rough sense, these remaining unresolved cases are not generic. 

As noted in the summary above, there are a number of other classes of free 
data for which we know whether or not solutions to (|47l) - (|48|l exist. Methods sim- 
ilar to those discussed here resolve the existence question for most CMC or near 
CMC free data which are asymptotically Euclidean or asymptotically hyperbolic, 
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as well for such data on closed manifolds. The situation for free data specified on 
compact manifolds with boundary is much less understood; however this is most 
likely a result of neglect rather than difficulty. Since numerical relativity has mo- 
tivated a recent interest in this question, and since the methods discussed here are 
believed to work for data on manifolds with boundary, it is likely that in the next 
few years, the solvability of (|47|l - 1)48(1 for such data is likely to be relatively well 
understood. Indeed, this motivation has led to the very recent results of Maxwell 
[?>T] and Dain j^H], who give sufficient conditions on asymptotically Euclidean free 
data with interior "apparent horizon" boundary conditions for solutions to exist 

The situation is very different for non-constant mean curvature data with no 
controls on |Vr|. Almost nothing is known, and there are no promising techniques 
known at this point. New ideas are needed. 

It is worth noting that the underlying approach of specifying the mean extrin- 
sic curvature r is geometrically natural, since the converse problem, of finding a 
hypersurface of prescribed mean curvature in a given spacetime, leads to a quasi- 
linear elliptic equation bearing some similarities with the Euclidean minimal surface 
equation. It is known that solutions of the Lorentz mean curvature equation are 
strictly spacelike and smooth, at least to the extent permitted by the regularity of 
the ambient spacetime and boundary conditions, so such prescribed mean curvature 
hypersurfaces provide natural spacelike slices of the spacetime [27| ^3 021 HOI ITT] . 
Although there are examples of spacetimes not admitting maximal (t = 0, |24|1 or 
CMC (r = const., HH) hypersurfaces, such non-existence behaviour is driven 
by global causal topology, which allows area-maximising sequences of hypersurfaces 
to become unbounded to the past or future; see |lllll2j . It is therefore not surprising 
that a priori control of the interior causal geometry is an essential ingredient in the 
proof of existence of entire maximal hypersurfaces in asymptotically flat spacetimes 

EH- 

To wind up this discussion of the implementation of the conformal method to 
find solutions of the constraint equations, it is important to note two facts: 1) The 
map from sets of free data to solutions of the constraints is surjective. This follows 
from the observation that if (7, K) is a solution of the constraints, then <j> = 1, W = 
is clearly a solution of (|47I48(I for the free data A = 7, a = K — ^Atr-yi^, r = tr 7 if . 
2) For all those sets of free data for which solutions of ((47148(1 are known to exist 
(excepting for the very special case of A flat, a = 0, and r = on M 3 = T 3 ) the 
solution is unique. 

As a consequence of these facts, those sets of free data for which solutions exist 
together parametrize the space of constraint-satisfying initial data for Einstein's 
equations. It follows that, once we determine exactly which free data sets map to 
solutions, we will have made significant progress toward understanding the "degrees 
of freedom" of Einstein's theory. 

Also as a consequence of these facts, if we seek initial data invariant under an 
isometry group, it is sufficient to choose free data with this invariance. 6 

In view of the dominant role the conformal method has assumed as a tool for 
mathematical as well as numerical analyses of the constraints, it is important to 
point out its limitations. As noted above, the conformal method replaces the original 
underdetermined PDE system of the constraint equations (four equations to be 
solved for twelve unknowns), by a determined system of elliptic character. This 
is very useful for a wide variety of studies (such as the parametrization question) . 
However for a number of other problems, it is a bad idea. Consider for example a 
given solution of the constraint equations on a finite radius ball. Can one smoothly 
extend the solution onto a neighborhood properly containing the ball? Can it be 

"Note that the presence of symmetries generally introduces integrability conditions which must 
be satisfied by the initial data, and consequently by the free data. These conditions are relatively 
easy to handle. 
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smoothly extended to a complete, asymptotically Euclidean solution? Might there 
be a smooth extension which is identically Schwarzschild or Kerr outside some 
(larger) ball? To study questions like these, casting the constraint equations into 
the form of a determined elliptic system is not at all useful. It is better to bypass 
the conformal method and work with the original underdetermined system. 



4.2 The Thin Sandwich Construction 

Some time ago, Wheeler El discussed a possible alternative to the initial value 
formulation for determining a spacetime solution of Einstein's equations. Rather 
than specifying initial data in the form of a Riemannian metric 7 and a symmetric 
tensor K which satisfy the Einstein constraint equations, he asked whether one 
might specify a pair of Riemannian metrics h\ and and seek to find a (unique?) 
spacetime (V 4 ,g) in which hi,h,2 arise as the induced metrics on disjoint Cauchy 
surfaces in (V 4 ,g), and which satisfies the Einstein equations i.e. {V 4 ,g) is the 
spacetime solution "connecting" h\ and hi. This is the sandwich conjecture. 

Considering the analogous question in electrodynamics, however, we find that 
the uniqueness part of the Maxwell version of the sandwich conjecture is false. Hence 
there is reason to believe that the uniqueness assertion in the sandwich conjecture 
for the Einstein equations may also be false. This led Wheeler to instead propose the 
thin sandwich conjecture, which postulates that given a freely chosen Riemannian 
metric 7 and symmetric tensor J representing the time derivative of the metric 
evolution in a spacetime, then we can use the constraints to determine a lapse 
function and shift vector field so that indeed, with respect to a choice of foliation 
and coordinates compatible with that lapse and shift, J is the time derivative of 
the metric. More explicitly, having chosen 7 and J, we seek a lapse N and shift X 
so that if we set 

K ij = -N-\y ij -X (i . j) ), (65) 

then (7, AT) satisfy the constraint equations (|19I20|) with source terms p — 167rT nn , 
Si = 8irT n i. Regarding (7, J, p, S) as given fields, the constraints give four equations 
for (N,X), the so-called thin sandwich equations. 

It is clear that if (7, K) satisfy the constraint equations, then J is the time 
derivative <9 t 7 of the spatial metric in the corresponding spacetime evolution of the 
data (7, if) with respect to coordinates with lapse-shift (N,X). 

The Hamiltonian constraint is readily solved for the lapse, giving 



where we have set T^- = \Jij — X(A-j\ for brevity. Substituting this expression for N 
into the momentum constraint l|2U|) we obtain the reduced thin sandwich equations 
("RTSE") 



V, 



(T*)2 - T kl T 



T5M 

(r«-r£ 7 «) 



= S\ (67) 



which we view as a system of partial differential equations for X 2 , the components 
of the shift. 

The linearization of (|67() is shown to be elliptic in [TJI, provided that the con- 
jugate momentum 7r y (recall the formula Q37(l for n) is either positive or negative 
definite. This rather surprising condition is guaranteed, so long as the condition 

p - R M > (68) 
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is satisfied everywhere in M and provided that J is chosen in such a way that the 
formula Ij66(l for the lapse function satisfies 

N > 0. (69) 

An additional necessary condition for solvability of the RTSE is that the equation 

= /xify (70) 

only has solutions {Vi,fi) with Vi — 0. 

An implicit function theorem argument |17| shows then that the RTSE is solv- 
able for all data (7, J, p, S) in an open neighbourhood of a reference configuration 
(7, J, p, S) satisfying the conditions (|68I69I70II . Here "open" is with respect to a 
Sobolev space in which the RTSE are posed. 

The fact that the thin sandwich equations can be solved for an open class of 
data (7, J, p, S) can be interpreted as showing that this reformulation of the con- 
straints leads to "generic" solutions. However, (|68|1 is rather restrictive — it excludes 
asymptotically flat data for example — so these conditions, taken together, limit 
the prospects of constructing most interesting spacetime initial data via the thin 
sandwich approach. 

4.3 Conformal Thin Sandwich 

An interesting hybrid of the conformal method and the thin sandwich approach has 
been suggested by York j76| . Before specifying how this conformal thin sandwich 
approach works, let us briefly compare the two approaches discussed so far: 

1. Conformal Method: The free data consists of a conformal metric Ay , a 
divergence- free trace-free tensor field cr^, and a scalar field r (8 free func- 
tions) . Solving the constraints produces a vector field W and a scalar field <j> 
(4 functions). We recompose to get a metric 7^ and a symmetric tensor Kij 
satisfying the constraints (12 functions). The lapse and shift are ignored. 

2. Thin Sandwich: The free data consists of a metric 7^ and a symmetric tensor 
Jij (12 free functions). Solving the constraints produces a vector field X 1 
and a scalar field N (4 functions). We recompose to get a metric 7^ and 
a symmetric tensor satisfying the constraints, plus the lapse N and the 
shift X 1 (16 functions). 

The idea of the conformal thin sandwich approach is to specify as free data 
a conformal metric Ajj, a trace- free symmetric tensor Uij, a scalar field r, and 
another scalar field r\. In terms of the loose function counting system used above, 
this amounts to 12 free functions. In the spirit of the thin sandwich, the tensor 
Uij represents the time derivative of the conformal metric, while r is the mean 
curvature and rj represents the lapse function, up to a conformal factor. It follows 
from this interpretation, together with a choice of conformal scaling, that the second 
fundamental form is expressed as 

Kio = ^±{{LX) l3 - Utj) + I^tXh (71) 

where X 1 is the shift vector and ip is the conformal factor, neither of which is known 
at this stage. To determine X % and ip, we use the constraint equations. These take 
the form 

V^((2r ? )- 1 (iA). u ) = V^r 1 ^-) + l^ViT, (72) 

and 

8AV> - R(-f)tp = - Aij A ij ip~ 7 + \t 2 ^. (73) 
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where we use the convenient short hand Aij — ^((LX)ij — Uij). Thus, if we can 
solve H72I730 f° r X and ip (4 functions), we have K (obtained from iJEJ) arid 

Hj = (74) 

satisfying the constraints, plus a specification of the lapse N = ip 6 rj and the shift 
X\ 

Compared to the conformal method ("CM"), the conformal thin sandwich ap- 
proach ("CTSA") has a number of virtues: 

1. Unlike CM, specifying a set of CTSA free data does not require a projection 
to the divergence- free part of a symmetric tracefree 2-tensor. 

2. The map from CTSA free data to a solution of the constraint equations is 
conformally covariant in the sense that if (7y, -Ky) ls a solution corresponding 
to a particular choice of free data (Ay, f7y,T, rf), then it is also a solution for 
the CTSA data (6 |4 Ay, 9~ 2 Uij, r, 6 r/), for any positive function 0. This is true 
whether or not r is constant, whereas CM data has this property only if t is 
constant. 

3. The mathematical form and hence the mathematical analysis of the CTSA 
equations f72 |) -l(75 )l is very similar to that of the CM equations |@7J| |@HJ|. 

Hence we have essentially the same existence and uniqueness results for the 
two sets of equations. 

4. CTSA free data is arguably closer to the physics we wish to model, since it 
includes the time derivative of the conformal metric, while the CM data only 
has the divergence-free trace-free (transverse traceless) part of the second 
fundamental form. 

Along with these virtues, CTSA has one troubling feature: Say we want to 
find a set of CMC initial data (7^, -Ky) with the lapse function chosen so that 
the evolving data continues to have constant mean curvature. In the case of the 
conformal method, after solving l|47|l and i|48|) to obtain a solution of the constraints, 
one achieves this by solving a linear elliptic equation for the lapse function. This 
equation is not coupled to the constraint equations, and solutions are readily verified 
to exist. By contrast, in the CTSA the extra equation takes the form 

A(V 7 ^) = t^V* + U^r\LX - Uf + ifx'ViT - V 5 , (75) 

which is coupled to the CTSA equations (|72Jl -173 (l . This coupling cannot be removed 
by choosing r constant. Indeed, whether the data is CMC or not, the mathematical 
analysis of the full system (|72H75|I is not very tractable. We emphasise, however, 
that if we do not require that the lapse be chosen to preserve the mean curvature, 
then the analysis of the CTSA equations is no more difficult than that of the con- 
formal method equations. For further discussion of the conformal thin sandwich 
approach, both in theory and in practice, see [TCj and |35| . 

4.4 Parabolic methods and the quasi-spherical ansatz 

It was first shown in |14| that 3-metrics of prescribed scalar curvature can be con- 
structed by solving a certain parabolic equation on S 2 . This leads to a method for 
constructing solutions of the Hamiltonian constraint which has found some inter- 
esting applications [1411711 R)§| . 

7 Although we restrict discussion here to 3-manifolds, the parabolic method generalizes to all 
dimensions. 
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To describe the general construction [71| . consider the second variation identity 
for a foliation r — ► S r in a Riemannian 3-manifold with metric g, 

R(g) = 2D n H + 2K G - H 2 - \\Ilf - 2u~ 1 Au , (76) 

where n is the outward unit normal vector to the level sets S r , II is the second 
fundamental form, H = trs r 77 is the mean curvature, and Kg, A = r~ 2 A are 
respectively the Gauss curvature and Laplacian of the metric r 2 g on the foliation 
2-surfaces (which are usually assumed to be topological 2-spheres). This is the 
Riemannian version of the Lorentzian formula i|35|) . The name "second variation" 
arises from the term D n H, since the mean curvature measures the first variation 
of the area of E r . This interpretation is not important for the present application, 
although it does reflect some deep relationship between mean curvature and minimal 
surfaces on the one hand, and positivity properties of mass on the other |^7| . 

Consider (|76() when the surface variation r — > E r defines a foliation, so the 
metric takes the general form 

g = u 2 dr 2 + g AB {rd9 A + (3 A dr){rd9 B + f3 B dr), (77) 

where u — u(r, 9 A ) > 0, (3 = {(3 A {r, 9 B )) and the rescaled angular metric gAB(r, 9 C ) 
are arbitrary fields. The quasi- spherical (QS) case considered in ^2] arises as the 
special case where g = dd 2 + sin 2 $ dip 2 is the standard 2-sphere metric. Generally, 
the mean curvature of the leaves S r is 

H = — ( 2 - div [3 + rd r log Jdet g) , (78) 
ru \ J 

where div /3 is the divergence of the tangential vector field (3 = (3 a 8a in the metric 
g, and n — u^ 1 (d r — r~ 1 f3 A dA) is the exterior-directed unit normal vector. Assume 
throughout the quasiconvexity condition 

h:=2- div o/ 3 + rd r log ^/det§ > 0, (79) 

and let n a — un, H a — uH = —h/r, II a = ull be the corresponding quantities 
determined by the metric H77JI with u = 1. The key observation is that substituting 
H from l|78|) into D n H gives a term hu~~ 3 r~ 1 d r u which combines with the term 
ti^'Au to show that if the scalar curvature is specified then l(7H|l may be read as a 
parabolic equation for the lapse u. Here (3 A , §ab can be freely chosen, subject only 
to the quasiconvexity condition 1)79(1 . or equivalently, H < 0. Explicitly we have 

h(rd r - f3 A d A )u - u 2 A D u (80) 
= ru(rd r - f3 A d A )H + \r 2 u{H 2 + \II \ 2 ) - u 3 (K(g) - ±R(g)r 2 ), 

which clearly shows the parabolic structure of the equation satisfied by u. 

Thus if R{g) is determined by the Hamiltonian constraint equation l|19f> and (3, 
c/ab are regarded as prescribed fields, then solving l(8()|l leads to a 3-metric which 
satisfies the Hamiltonian constraint. 

Global (r — > oo) existence theorems were established for l|80|l in the quasi- 
spherical case g = dd 2 + sin 2 •& dip 2 ^1] , assuming (3 satisfies the quasiconvexity 
condition (|79|l . which becomes div G /3 < 2, and assuming the prescribed scalar cur- 
vature is not too positive (r 2 R < 2 is sufficient but not necessary). With suitable 
decay conditions on the prescribable fields (3 A (for example, if (3 A is compactly 
supported), these results give asymptotically flat solutions of the Hamiltonian con- 
straint, with either black hole H = or regular centre r — inner boundary condi- 
tions. Setting R(g) = gives a large class of solutions of the vacuum Hamiltonian 
constraint l(19|l in the time- symmetric case — 0, with two degrees of freedom 
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corresponding to the choice of (3 A . Analogous results have been established for 
general Ijab |71l 1721 which demonstrate that this provides a flexible technique 
for extending solutions to be asymptotically flat, or asymptotically hyperboloidal if 
R(g) = -6. 

The parabolic method enables us to construct metrics of prescribed scalar cur- 
vature which have properties not achievable by conformal methods. For example, 
any bounded domain (fl, g) with smooth mean-convex (H < 0) boundary dfl ~ S 2 
and non-negative scalar curvature can be extended smoothly to an asymptotically 
flat manifold, also with R > 0; we simply extend the Gaussian normal foliation in 
a neighbourhood of dQ smoothly to a metric which is flat outside a neighbourhood 
of Q, with a mean-convex exterior foliation approaching the standard spherical foli- 
ation of R 3 . Solving the parabolic lapse equation for R > with initial data u — 1 
on a Gaussian level set £_ e C ft produces the required 3-metric. 

One consequence of this argument, first observed in |14j . is the existence of 
R = metrics on M 3 which are asymptotically flat and non-flat, but which have a 
flat interior region. 

The parabolic method also gives solutions satisfying "geometric boundary condi- 
tions" |16j (r 2 <7, H), which specify the boundary metric t 2 c/ab and mean curvature 
H. This is clear from l|78|l since H is determined by specifying the well-posed ini- 
tial condition u(ro) for the parabolic lapse equation. Note that this is not possible 
with conformal methods, since specifying both the boundary metric and mean cur- 
vature leads to simultaneous Dirichlet and Neumann boundary conditions on the 
conformal factor. These boundary conditions are ill-posed for the elliptic equation 
imposed on the conformal factor, which will not then admit any solution in general. 
The problem of constructing 3-metrics with prescribed boundary metric and mean 
curvature arises naturally from the geometric definition of quasi-local mass |13| and 
provided the original motivation for the quasi-spherical construction. 

Examples of 3-metrics satisfying the Geroch (inverse mean curvature flow) foli- 
ation condition ^| QH| can be constructed by choosing (3 A , §ab satisfying div D /3 = 
rd r (log det g) and then solving the parabolic lapse equation with suitably prescribed 
scalar curvature. This provides a large class of metrics for which the Geroch iden- 
tity directly verifies the Penrose conjecture. Note that this does not help to solve 
the much harder problem considered in |48| . of finding a Geroch foliation in a given 
metric. 

Another advantage of the parabolic lapse method is that it is generally easier 
and cheaper numerically to solve a 2+1 parabolic equation than the 3D elliptic 
equations arising in the conformal methods, particularly when the solution is re- 
quired near spatial infinity. Finally, in the original quasi-spherical gauge |14j . the 
freely-specifiable fields (3 A ,A = 1,2 give a rather explicit parameterization of the 
"true degrees of freedom" allowed by the Hamiltonian constraint, since fixing j3 
removes all diffeomorphism freedom in the quasi-spherical gauge (at the linearized 
level at least). 

Dual to the problem of constructing metrics having quasi-spherical form, is 
the harder problem of finding a quasi-spherical foliation in a given metric. Some 
simple observations suggest that QS foliations exist for general perturbations of a 
QS metric, but for interesting technical reasons a complete proof is not yet available. 

Write the general QS 3-metric as g — u 2 dr 2 + E^/^dr + toa) 2 where <ja, 
A = 1, 2 is an orthonormal coframe for the unit sphere S 2 , and consider the effect 
of an infinitesimal diffeomorphism generated by X — C ta + z(rd r — (3 a ta), where 
ta is the dual frame to a a- The QS condition will be preserved under the metric 
variation Sgij and infinitesimal diffeomorphism X exactly when (Cx9 + $9)ab = 0, 
which gives the equation 

C(A\B) + {°AB ~ P(A\B))z + \5g A B = 0. (81) 



24 



Taking trace and trace-free projections and simplifying gives 

VpC := C(A|B) - ^div (<r AB + B AB div ( 

= -\{&9ab - \&9 C C <JAB + B AB 6g%), (82) 

where 

„ P(A\B) ~ \d\v t3a A B 
2 - div /3 

and 

= div C + \5gj 
2-div D /3 

The operator VpC, in (|82|) is elliptic exactly when \B\ 2 < 1, and we can show that 

its adjoint T>ly AB = —y AB + {ubcB bc )\ A has trivial kernel when \B\ 2 < 1/3 
(pointwise). This latter condition also guarantees that Dp is surjective, whereupon 
(|82|l can be solved for ( A , uniquely if the L = 1 spherical harmonic components 
of xa are specified. Hence existence of QS coordinates can be established at the 
linearized level, but the presence of div Q £ in the expression for z means that it is not 
possible to directly apply an implicit function theorem to conclude the existence for 
all metrics in an open neighbourhood of g. 

It is also interesting to consider the direct approach to finding quasi-spheres as 
graphs over S 2 in M 3 . From the second fundamental form of the graph r = e u , 
u : S 2 -> M, 

e u . , 

II AB = , j =\ U AB - U A U B - <J AB ) 

y/1 + \DU\ 2 

where a AB is the S 2 metric and |Z?u| 2 = cr AB u A u B , the Gauss and mean curvatures 
are given by 

K(u) = det II AB j det a AB (83) 
H{u) = e- 2u a AB II AB , (84) 

where a AB = e 2u (a AB + 11,4115) is the induced metric on the graph. It follows 
that the linearisations about the unit sphere with respect to infinitesimal changes 
v = Su, h = Sg are 

-5K(v,h) = (A + 2)v + ±((A + l)(a AB h AB )-h AB lAB ) (85) 
Hv = (A + 2)v + h rr + a AB (h AB + S7 A h Br -^d r h AB ). (86) 

Now the equation (A + 2)<fi — f is solvable for 4> exactly when <f g2 fY\ = for any 
I = 1 spherical harmonic Y\, and it follows that for generic metric variations h, it 
is not possible to find a corresponding variation v of the graph which preserves the 
condition H = — 2. This instability of constant mean curvature foliations can be 
overcome if the background metric has positive ADM mass I n contrast, the 

equation SK(v, h) = is solvable for v for any metric variation h, since 

/ (A + l)a AB h AB - h l ^)Y 1 = 

for all Y\ satisfying AYi = — 2Y\ and then VYi is a conformal Killing vector of S 2 
and V^ B Yi = (AYi + Yi)a AB . Thus at the linearised level it is always possible to 
preserve the condition K = 1. 

Further evidence supporting the conjecture that QS foliations exist for generic 
metrics near a given QS metric satisfying \B\ 2 < 1/3 comes from the thesis of 
Spillane |73j . which establishes existence in the axially symmetric case. However, 
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the general case remains open, as does the related problem of finding a QS foliation 
of a null hypcrsurface ^Hj- Note that this foliation existence question is irrelevant 
for the metric existence question, which is resolved by the proof in |14j of a large 
class of metrics in QS form. 

The momentum constraint equations are not yet well-understood in the parabolic 
method, although some results have been established 6&. In the quasi-spherical 
case, we may introduce the parameterization 

K ij i 6 j = (r] AB + ±(iS AB )e A 8 B + n A {8 A 8 3 + 6 3 6 A ) + (r - f i)6 3 8 3 , (87) 

where 6 1 = (3 1 dr + rd/&, 9 2 = (3 2 dr + rs'm-ddip, 9 3 = udr is the QS orthonor- 
mal coframe. In terms of the parameters t]ab, (J-, k-a, t = tr^K, the momentum 
constraint equations then take the form 

87rT 03 ™ = -(rd r ju-/3 A VA/«+(3-div/3)/i) (88) 
+ udiv k + k a X7 a u + T] AB 0(a\b) + (2 - div 0)t 

8nT 0A ru = (ut] A b) ]B + ^Va« + Va(u(|p - T )) (89) 
+ rd r KA - P B ka\b + ((3 - div 0)8 ab ~ Pb\a)k b ■ 

The matter fields Too, Tqa, T33 are prescribed and l|88|l provides either an equation 
for r, by choosing /1, or vice-versa, choosing r and solving for fi. Likewise, (I89f) can 
be regarded either as an equation for n (with the symmetric traceless 2-tensor tjab 
arbitrarily prescribable) or as an elliptic equation for rj with ka freely prescribable. 

Local existence in r for the momentum constraints with prescribed r, ka has 
been established by Sharpies (HHjj but it is not clear whether global results are 
possible without additional restrictions. Much work still remains to be done on 
these systems. 

Finally we note that a characteristic version of the quasi-spherical gauge, where 
a foliation of outgoing {pnp < 0) null hypersurfaces is assumed to admit a QS radial 
coordinate, has been described in |15| . In this case there is no parabolic equation. 
The resulting hypersurface Einstein equations are considerably simpler than those 
derived in the Bondi gauge 66, 23 , and forms the basis for a 4th order numerical 
code ^Sjj which heavily exploits the exact spherical geometry of the r-level surfaces. 

4.5 Gluing Solutions of the Constraint Equations 

The conformal method, the conformal thin sandwich method, and the quasi-spherical 
ansatz are all procedures for generating solutions of the Einstein constraint equa- 
tions from scratch. We now consider procedures for constructing new solutions of 
the constraints from existing solutions. 

We first discuss a procedure for gluing connected sums of solutions. The idea of 
this "IMP Gluing" is the following: Say we have two solutions of the constraint 
equations, (Mi, ji,Ki) and (M2, 72, -^2)- Let pi £ M% and p 2 S M 2 . Can we 
find a set of initial data (Mn_ 2 ), 7(i-2)> -KV1-2)) such that 1) Mi_ 2 is homotopic 
to the connected sum 8 Afi#M 2 ; 2) (7(1-2)5-^(1-2)) is a solution of the constraints 
everywhere on Mn_2); and 3) on that portion of Mn_2) which corresponds to 
Mi \ {ball around pi}, the data (7(1-2)1^(1-2)) is isomorphic to (71,^1), with a 
corresponding property holding on that portion of M(i_ 2 ) which corresponds to 
M 2 \ {ball around p 2 }? If so, we say that the sets of data admit IMP gluing. 

IMP gluing can be carried out for quite general sets of initial data. The sets 
can be asymptotically Euclidean, asymptotically hyperbolic, specified on a closed 

8 The connected sum of these two manifolds is constructed as follows: First we remove a ball 
from each of the manifolds Mi and M2. We then use a cylindrical bridge S 2 X I (where / is an 
interval in R 1 ) to connect the resulting S 2 boundaries on each manifold. 
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manifold, or indeed anything else. The only condition the data sets must satisfy is 
that, in sufficiently small neighborhoods of each of the points at which the gluing is 
to be done, there do not exist nontrivial solutions £ to the equation D$*^ K ^ = 0, 
where D<&J K ^ is the linearization operator defined in i|43|) (with K replacing 7r). 
In |2(J| it is shown that this condition is indeed generic. 

The proof that IMP gluing can be carried out to this degree of generality is 
detailed in |33|, based on j52], jSH], and We note here two features of it. 

First, the proof is constructive, in the sense that it outlines a systematic, step- 
by-step mathematical procedure for doing the gluing: One conformally blows up 
the balls surrounding p\ and pi to produce two half cylinders extending from the 
original initial data sets; one joins the two half cylinders into a bridge, and splices 
together the data from each side using cutoff functions; one uses the local constant 
mean curvature to decouple the constraints in the neighborhood of the bridge; 
one uses tensor projection operators based on linear PDE solutions to find a new 
conformal K which solves the momentum constraint; one solves the Lichnerowicz 
equation (the argument that this can be done, and that the solution is very close to 
1 away from the bridge, relies on the invertiblity of the linearized equation and on 
a contraction mapping); one recomposes the data as in H50|l and (|51(l : and finally 
one does a nonconformal data perturbation away from the bridge to return the 
data there to what it was before the gluing. This procedure can be largely carried 
out numerically, although we note that it requires us to solve elliptic equations on 
topologically nontrivial manifolds. 

The second feature we note regarding the proof is that it relies primarily on the 
conformal method, but it also uses a nonconformal deformation of the data at the 
end, to guarantee that the glued data is not just very close to the given data on 
regions away from the bridge, but is exactly equal to it. 

While IMP gluing is not the most efficient tool for studying the complete set of 
solutions of the constraints, it has already proven to be very useful for a number of 
applications, including the following: 

1. Multi-Black Hole Data Sets: Given an asymptotically Euclidean solution of 
the constraints, IMP gluing allows a sequence of flat space initial data sets to 
be glued to it. The bridges that result from this gluing each contain a minimal 
surface, and consequently an apparent horizon. With a bit of care |32| . one 
can do this in such a way that indeed the apparent horizons are disjoint, and 
therefore likely to lead to independent black holes. 

2. Adding a Black Hole to a Cosmological Spacetime: Although there is no clear 
established definition for a black hole in a spatially compact solution of Ein- 
stein's equations, one can glue an asymptotically Euclidean solution of the 
constraints to a solution on a compact manifold, in such a way that there is 
an apparent horizon on the bridge. Studying the nature of these solutions of 
the constraints, and their evolution, could be useful in trying to understand 
what one might mean by a black hole in a cosmological spacetime. 

3. Adding a Wormhole to Your Spacetime: While we have discussed IMP gluing 
as a procedure which builds solutions of the constraints with a bridge connect- 
ing two points on different manifolds, it can also be used to build a solution 
with a bridge connecting a pair of points on the same manifold. This allows 
one to do the following: If one has a globally hyperbolic spacetime solution 
of Einstein's equations, one can choose a Cauchy surface for that solution, 
choose a pair of points on that Cauchy surface, and glue the solution to itself 
via a bridge from one of these points to the other. If one now evolves this 
glued-together initial data into a spacetime, it will likely become singular very 
quickly because of the collapse of the bridge. Until the singularity develops, 
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however, the solution is essentially as it was before the gluing, with the ad- 
dition of an effective wormhole. Hence, this procedure can be used to glue a 
wormhole onto a generic spacetime solution. 

4. Removing Topological Obstructions for Constraint Solutions: We know that 
every closed three dimensional manifold M 3 admits a solution of the vacuum 
constraint equations. To show this, we use the fact that M 3 always admits 
a metric T of constant negative scalar curvature. One easily verifies that the 
data (7 = r, K = T) is a CMC solution. Combining this result with IMP 
gluing, one can show that for every closed M 3 , the manifold M 3 \ {p} admits 
both an asymptotically Euclidean and an asymptotically hyperbolic solution 
of the vacuum constraint equations. 

5. Proving the Existence of Vacuum Solutions on Closed Manifolds with No CMC 
Cauchy Surface: Based on the work of Bartnik |1 11 112) . one can show that if 
one has a set of initial data on the manifold T 3 #T 3 with the metric compo- 
nents even across a central sphere and the components of K odd across that 
same central sphere, then the spacetime development of that data does not 
admit a CMC Cauchy surface. Using IMP gluing, one can show that indeed 
initial data sets of this sort exist. 



4.6 The Corvino-Schoen method 

There is another very useful form of gluing which has been applied recently to con- 
struct interesting solutions of the Einstein vacuum constraint equations. Developed 
by Corvino and Schoen [361 1371 l3~2~] , this method has the following remarkable ap- 
plication. Let (M 3 ,7,-ftT) be a smooth, asymptotically Euclidean solution of the 
constraint equations. If certain asymptotic conditions hold, then for any compact 
region E 3 C M 3 for which M 3 \ S 3 = M 3 \ S 3 (where B 3 is a ball in M 3 ), there is a 
smooth asymptotically Euclidean solution on M 3 which is identical to the original 
solution on T, 3 C M 3 , and is identical to Cauchy data for the Kerr solution on 
M 3 \ S 3 for some S 3 C M 3 . In words, their technique allows us to smoothly glue 
any interior region of an asymptotically Euclidean solution to an exterior region of a 
slice of the Kerr solution. For asymptotically Euclidean solutions of the constraints 
with tXjK — 0, this method glues any interior region to an exterior region of a slice 
of Schwarzschild. 

Combining the Corvino-Schoen gluing techniques with some results of Friedrich 
44 , 32 showed that there is a large class of vacuum spacetime solutions of Ein- 
stein's equations which admit complete null infinity regions of the form "scri" , as 
hypothesized by Penrose. The tools developed by Corvino and Schoen have also 
been used to strengthen the IMP gluing results |33| . 

The Corvino-Schoen method aims to solve the constraint equations through a 
projection using the linearized operator _D$ and its adjoint D$>*. We sketch the 
method in the time-symmetric case 7r = 0, where $(7, 7r) is replaced by the scalar 
curvature R(j) and the lapse-shift £ is replaced by the lapse N. In this case the 
arguments are essentially the same while the calculations are considerably simpler. 

We start with the observation that because DR* has injective symbol, it satisfies 
an elliptic estimate on any domain 0, 

\\N\\ HHn) < C(\\DR*N\\ L 2 (n) + \\N\\ mn) ), (90) 

which importantly does not require any control on N at the boundary dQ. It follows 
easily that similar weighted estimates hold, with weight function p £ C^° which is 
positive in f2 and vanishes to high order at dfl: 

\\N\\m<n) := / P(N 2 + |V 2 iV| 2 ) dv, < C f p\DR*N\ 2 + N 2 ) dv„ (91) 
" Jn Jn 
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where the final term N 2 on the right can be removed if there are no Killing vectors. 
With this assumption, for all / G L 2 (Ct) we can solve DR(pDR*N) = pf for 
N G Hioc wmcn m particular produces a solution to the linearized constraint 
equation DRh = pf. An iteration argument is used to solve the nonlinear problem 
R(jo + h) — R(jo) + S for any sufficiently small S. This solution h G H 2 _ 1 has 
the remarkable property that it vanishes to high order on dil. Thus, for example, 
if i?(7o) is sufficiently small and supported in Q then there is a perturbation h, also 
supported in f2, such that R("fo + h) =0. 

To use this method to glue a Schwarzschild exterior to an asymptotically flat 
R(~f) — metric across an annulus B2r\Br, R >> 1, requires one more idea 
because the flat space kernel keiDRg = span(l, x , x 2 , x 3 ) is non-trivial. This 
implies that the linearized problem DR$h = a is solvable if and only if a satisfies 
the four conditions J n er(l, x l )d 3 x = 0, and the nonlinear problem R(jo + h) = is 
similarly obstructed for 70 close to flat. By choosing the cutoff radius R sufficiently 
large and rescaling back to ft = B2\B\ produces exactly this close to flat situation. 
However, it is possible to solve the projected problem i?(7o + h) G K := span(l, x % ) 
with uniform estimates on ft £ H 2 _ 1 {^1). Now the Schwarzschild exterior metric 
can be characterised by the mass and centre of mass parameters (to, c*), defined by 

7schw = (l+2^)%- (92) 

Some delicate estimates show that the map (to, c l ) 1— ► K is continuous and has 
index 1, so there is a choice of parameters (to, c l ) mapping to G K, which gives 
R(-fo + h) — as required. The extension of these ideas, and the considerable details 
of the above arguments, are given in the original references (SElEZIESj- 

5 Conclusion 

A considerable amount is known concerning the solutions of the Einstein constraint 
equations. Using the conformal method or the conformal thin sandwich method, we 
know how to construct constant mean curvature solutions which are asymptotically 
Euclidean, asymptotically hyperbolic, or live on a closed manifold. We also know 
how to do the same for nearly constant mean curvature solutions. We can glue 
together quite general solutions of the constraints, producing new solutions of both 
mathematical and physical interest. And, for certain asymptotically Euclidean so- 
lutions, we know how to show that there are solutions which include any compact 
region of the solution in the interior, and which are exactly Kerr or Schwarzschild 
in the exterior. 

Much remains to discover as well. We would like to know how to construct 
solutions with mean curvature neither constant nor nearly constant. We would like 
to know much more about constructing solutions of the constraints with prescribed 
boundary conditions. We would like to know to what extent we can construct solu- 
tions with low regularity. And we would like to know which solutions on compact 
regions can be smoothly extended to either asymptotically Euclidean or asymptot- 
ically hyperbolic solutions. 

Besides these mathematical issues to resolve, there are important questions con- 
cerning solutions of the constraints and physical modeling. In view of the pressing 
need to model astrophysical events which produce detectable amounts of gravita- 
tional radiation, one of the crucial questions we need to answer is how to system- 
atically find solutions of the constraint equations which serve as physically realistic 
initial data sets for such astrophysical models. Since these models are generally 
constructed numerically, an equally crucial question is the extent to which the con- 
straint functions, initially zero, remain near zero as the spacetimc is numerically 
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evolved. 

It is not clear how close we are to resolving these mathematical and physical 
questions regarding the Einstein constraint equations and their solutions. However, 
in view of the recent rapid progress that has been made in these studies, we are 
optimistic that many of them will be resolved soon as well. 
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